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Abstract 

In case of chiral Uyi(l) symmetry restoration the mass of the rj' boson (the 
ninth, would-be Goldstone boson) is decreased, thus its production cross 
section is heavily enhanced [1]. The rj' decays (through one of its decay 
channels) into five pions. These pions will not be correlated in terms of Bose- 
Einsten correlations, thus the production enhancement changes the strength 
of two-pion correlation functions at low momentum [2]. Preliminary results 
strongly support the mass decrease of the rj' boson [3H5]. In this paper we 
propose a method to select pions coming from rj' decays. We investigate the 
efficiency of the proposed kinematical cut in several collision systems and 
energies with several simulators. We prove that our method can be used in 
all investigeted collision systems. 



1 Chiral symmetry and rj' mass 

The temperature of the quark gluon plasma created in gold-gold collisions of 
the Relativistic Heavy Ion Collider (RHIC) may reach values up to 300-600 
MeV [6113 ■ these very high temperatures the degrees of freedom are not 
hadrons but quarks or gluons. It is expected that the broken symmetries 
of QCD are partially restored in this matter. Originally, the Ua{^) part of 
chiral symmetry is broken exactly, thus a high mass meson is produced, the 
rj' particle, which has a mass of 958 MeV. However, if the Ua{^) symmetry 
is partially restored, the mass of the rj' is decreasing [1]. 

Production cross sections of hadrons are exponentially suppressed by their 
mass. Hence without mass modification roughly two orders of magnitute less 
rj' mesons are produced than pions. In contrast, decreased mass rj' mesons 
will be created more abundantly. Thus the number of rj' mesons is closely 
related to their mass. 
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The decay of the 77' happens after it regained its vacuum mass (at the 
expense of its momentum). One important decay channel is the decay into 
two leptons, 77' — )■ + this is investigated in ref. j8]. It turns out that 
there is an excess in the dilepton spectrum at low invariant mass, and this 
excess might be related to the 77' enhancement. There is also a decay mode 
when the v( goes into an 77 and two pions, and the i] also decays into three 
pions: 

77' — 7/ + 7r+ + tt" — (vr^ + tt" + 7r°) + 7r+ + tt^ (1) 

and the overall probability of this decay chain is 10% [9]. The average mo- 
mentum of the resulting five pions is 138 MeV due to the low momentum of 
the original r( [2]. 

2 Two-pion Bose-Einstein correlations 

Final state effects distort two-particle correlation functions. One of the most 
important final state effect is that of Bose-Einstein correlations. To calculate 
these correlation functions, let us utilize the core-halo model. In the core- 
halo model [TU], the hadronic source is divided into two parts: a core and 
a halo. The core consists of the primordial particles and decay products of 
very short lifetime resonances, thus its size is very small, roughly 5 fm. The 
halo consists of decay products of long lived resonances, such as 77, j]' or 
K^g. The halo hadrons are created very far from the core, so the halo size is 
much larger, at least 50 fm. When measuring correlation functions however, 
due to finite momentum resolution of the detectors very small momentum 
differences cannot be resolved, i.e. pairs with such similar momenta are 
regarded as one by the detectors. In the Fourier transformation, large sizes 
correspond to small momenta, thus the halo correlations are not seen in 
measurements. Hence the observable part of the correlation function is due 
to the core, but its strength is decreased by the core/halo ratio. In case of a 
plain-wave approximation, for an identical boson pair with momenta and 
P2, the two-particle correlation function is (for details see e.g. ref. |10)): 

C.(g,i^) = l + A >^^'^^'' , (2) 
|^c(g = o,ir)| 

where q = Pi —p2, K = {pi +p2)/2 and Sc{q,p) is the Fourier transformed of 
the core source function (or emission function of pions coming from the core) 
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Sc{x, K) (the Fourier transformation is in x — g). The important parameter 
here is A = Nc/{Nc + Nh) , where the number of particles in the core is Nq, 
number of particles in the halo is A'^^- The A parameter thus depends on the 
ratio of the core to the halo. The halo pions come partly from t]' decays, so 
more rj' means larger halo, and smaller A. Hence t]' mass and A value are 
connected. This was investigated in refs. [3H5], and it was found that the A 
parameter is indeed decreasing at the kinematical domain of rj' decay pions. 
However, it is not experimentally proven that the rj' decay pions are causing 
the decrease. In this paper we investigate a method to kinematically filter 
out pions from rj' decays. If applied to the experimental sample, in case of 
an rj' mass modification the A decrese will vanish. 

3 Kinematical domain of pions from rj' decays 

Our method is based on the invariant mass of pions from the decay chain 
mentioned in eq. [TJ Let us investigate the invariant momentum of tt'^,tt~ 
pairs in the above decays. Using the mass-shell condition for pions, and 
utilize momentum conservation in the decay, one gets an rriinv interval for 
the 7r~^,7T^ pair coming from the 77', for the second pion pair coming from 
the f], and for the whole quadruplet. Based on the simulations, we chose the 
0.075-0.171 GeVVc^ interval for pairs and the 0.43-0.69 GeVVc^ interval for 
quadruplets. These yields an effective method of kinematical selection of rj' 
decay pions, as detailed in the next section. 

Our method is the following. For a given 7r~^, we take all possible pion 
quadruplets (7r"^,7r~,7r"'',7r~) from the same event, and check if any of these 
quadruplets fulfill our rriiny criteria. If there is at least one, we consider the 
original vr"*" to be "found", or say that if fulfills our rriinv criteria. The same 
can be done by starting from a given 7r"^,7r~ pair. In this case "finding" a 
non-?7' pion pair will be less probable, but the efficiency of finding rj' pions 
will not be different. So this is expected to be a better working method. In 
simulations, we can determine if the pair of the particle comes from an rj' 
decay, so the efficiency of the selection method can be tested. Whether using 
the pair or the single particle method, we can form four different groups of 
them: 

a) Comes from an 77' and fulfills the rriinv criteria 

b) Comes from an 77' and does not fulfill the rriinv criteria 

c) Does not come from an 77' and fulfills the rriinv criteria 

d) Does not come from an 77' and does not fulfill the rriinv criteria 
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Figure 1: Results from ^/s = 200 GeV p+p collisions. The method is working 
in all cases, but efficiency is very low in case of PHENIX cuts, for the single 
particle method. 

Ideally all r]' pions go into the ffist group, while all others fall into the last 
group. Let us call the number of pions (or pairs in the other method) in 
the four groups Na, N^, N^. and N^, respectively. In the next part, we will 
give the following numbers for each type of simulation: Efficiency of cutting 
rj' decay products, Na/ {Na + Nb) ; and loss of statistics, i.e. fraction of lost 
non-?7' pions, N^/ {N^ + N^) ■ The optimal value for efficiency is 1 (in this 
case we could cut out all pions coming from 7]' mesons), for loss (in this 
case we kept all non-r/' pions). Note that a "found" pair or particle that is not 
coming from an rj' (i.e. a loss greater than zero) decreases our experimental 
sample. If the loss is 90%, then the sample is reduced by a factor of 10, so the 
statistical errors will heavily increase. Goal of present paper is to investigate 
the efficiency and loss connected to our the method. A similar method was 
investigated in ref. for e'^e~ collisions. We test the method in p+p and 
Au+Au collisions, at several center-of-mass energies. 
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Figure 2: Results from ^/s = 200 GeV p+p collisions. The method is working 
in all cases. 
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Figure 3: Results from y/s = 200 GeV Au+Au collisions. The method is 
working only in the pair method, as loss is 100% in the other case. 
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4 Results 



We used two simulations to test our method: Pythia [l^ (proton-proton col- 
lisions, version 8.135) and Hijing [13] (proton-proton and gold-gold collisions, 
version 1.411). We also simulated the geometric acceptance of the detectors. 
In case of the 200 GeV RHIC energy, we used the geometry of STAR and 
PHENIX detectors, while in case of 14 TeV energy, we used the geometry of 
ALICE and CMS detectors. 

We generated 1 000 000 ^/s = 200 GeV and 10 000 = 14 TeV p+p 
events. For PHENIX cuts, the efficiency is the lowest, because a large fraction 
of pions are not detected. At both energies, the pair method is better than 
the single particle method. We also generated 100 ^/sp^ = 200 GeV events. 
Here only the pair method is working, as essentially all single particles are 
found, due to the very large statistics. See details on figs. [T]l3] for details. 
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